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Abstract 

In this paper, which is part II in a series of two, the pre-asymptotic error analysis of the continuous 
interior penalty Hnite element method (CIP-FEM) and the FEM for the Helmholtz equation in 
two and three dimensions is continued. While part I contained results on the linear CIP-FEM 
and FEM, the present part deals with approximation spaces of order p > 1. By using a modified 
duality argument, pre-asymptotic error estimates are derived for both methods under the condition of 

< Co{^) ''^^ , where k is the wave number, h is the mesh size, and Co is a constant independent of 
k, h,p, and the penalty parameters. It is shown that the pollution errors of both methods in H^-norm 
are 0{k'^P+^ h'^") if p = 0(1) and are ©(^(l^)^") if the exact solution u € H^{n) which coincide 
with existent dispersion analyses for the FEM on Cartesian grids. Here cr is a constant independent 
of k,h,p, and the penalty parameters. Moreover, it is proved that the CIP-FEM is stable for any 
k,h,p > and penalty parameters with positive imaginary parts. Besides the advantage of the 
absolute stability of the CIP-FEM compared to the FEM, the penalty parameters may be tuned to 
reduce the pollution effects. 

Key words. Helmholtz equation, large wave number, pre-asymptotic error estimates, continuous 
interior penalty finite element methods, finite element methods 
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1 Introduction 

This is the second installment in a series (cf. [17]) which is devoted to pre-asymptotic stability and error 
estimates of hp versions of some continuous interior penalty finite element method (CIP-FEM) and the 
finite element method (FEM) for the following Helmholtz problem: 

(1.1) -Au- k^u = f infl, 

Ou 

(1.2) — + iku^ g on F, 

on 

where C K'^, d = 2, 3 is a domain with smooth boundary, F :— 9f2, i = denotes the imaginary 

unit, and n denotes the unit outward normal to dQ. The above Helmholtz problem is an approximation 
of the following acoustic scattering problem (with time dependence e"^*): 

(1.3) -Au-k^u^f inM'*, 

(1.4) ( ^^"^ + ik{u - m'"")) ^ as r = |a;| ^ C50, 

where u™'^ is the incident wave and k is known as the wave number. The Robin boundary condition 
(|1.2p is known as the first order approximation of the radiation condition (|1.4[) (cf. [H]). We remark 
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that the Helmholtz problem (|l.ip - (|1.2[) also arises in applications as a consequence of frequency domain 
treatment of attenuated scalar waves (cf. p3)). 

It is well-known that the error bound of the finite element solution to the Helmholtz problem (jl.ip - 
(|1.2[) usually consists of two parts, one is the same order as the error of the best approximation of u from 
the finite element space, another one is worse than the convergence order of the best approximation and 
dominates the error bound of the finite element solution for large wave number k (relative to h and p) 
[D[5l[i[2TJ[2g[2i[23|3Tl[3l|33|45l|47]. The second part is the so-called poUution error in the literature. 
We recall that, the term "asymptotic error estimate" refers to the error estimate without pollution error 
and the term "pre-asymptotic error estimate" refers to the estimate with non-negligible pollution effect. 
The latest asymptotic error analysis was given by Melenk and Sauter [39l|40]. It is shown that, the 
hp-FFiM is pollution-free in the i7^-norm if 

(1.5) 1- k(—Y is small enough, 

p ap 

where h is the mesh size, p is the polynomial degree of finite element space, and a is some positive constant 
independent of k, h, and p. In particular, the hp-FFM is pollution-free under either of the following two 
conditions: 

(1.6) p — 0(1) fixed independent of k and k^^'p h is small enough. 



(1.7) p> clnk and — is small enough, 

P 

where c is some positive constant independent of fc, h, and p. Although the above results improve greatly 
the previous results which require small enough (cf. [331 [3Sj), the pre-asymptotic error estimates are 
still worth to be studied based on but not limited to the following considerations: 

• For fixed p, the pollution effect can be reduced substantially but cannot be avoided in principle (cf. 

mm)- 

• Implementations of high order FEMs, e.g. p > Ink for large fc, are not easy for problems with 
complicated geometries or different materials. 

• Given a tolerance e, one wish to use as small as possible number of degrees of freedom, or in other 
words, as large as possible mesh size h, to achieve this tolerance. We denote the maximum mesh 
size by h{k,e). Then h{k,e) always locates in the pre-asymptotic range for large k because the 
pollution error dominates the error bound. For example, when p — 1, for large k, the pollution 
term is 0{k^h'^) (cf. [331 SZ]) and hence h{k,e) — ©((e/fc^)^/^) which is not in the asymptotic 
range (0,C/fc^) implied by the condition (|1.6p . 

• The pre-asymptotic error analysis of FEM is hard for higher dimensional problems. To the best 
of the authors' knowledge, besides the first part of this series [37] for the linear FEM, no work 
has been done ever since the pioneer works of Ihlenburg and Babuska [34[ 135] for one dimensional 
problems. 

The purpose of this paper and the companion paper |47j is twofold: One is to derive pre-asymptotic 
error estimates for the FEM for the Helmholtz equation in two and three dimensions, under the condition 

of -y^ < C'o(f)''^^ for some positive constant Co independent of fc, ft,, and p. Clearly, this condition 
extends the asymptotic range given by ()1.5p . For example, when p = 0(1), the above condition is 
satisfied if /c^^p+r/i is small enough. This would be the condition obtained for the {p -\- l)-th order 
FEM by "standard" arguments (cf. (II. 6^ ). Another purpose is to analyze some CIP-FEM which is 
absolutely stable (that is, stable for any fc, /i, and p) and capable of achieving much less pollution effect 
than the FEM for fixed p. The CIP-FEM, which was first proposed by Douglas and Dupont [55] for 
elliptic and parabolic problems in 1970's and then successfully applied to convection-dominated problems 
as a stabilization technique [T31 [T31 UHl HH HZ] , modifies the sesquilinear form of the FEM by adding 
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the following least squares term penalizing the jump of the gradient of the discrete solution at mesh 
interfaces: 



(1.8) J(u,«):^^i7| 



du 
drip. 



dv 



Note that the penalty parameter 17 is chosen as a complex number in this paper and 47 instead of the 
real numbers as usual. The paper [47' has considered the linear case p = 1 and this paper will be devoted 
to the /ip-version. To be precise, we obtain the following results: 

(i) There exists a constant Co > independent of k, h, and 7, such that if fc > 1, < 7 < 1, and 



p+i 



then the following pre-asymptotic error estimates hold for both the CIP-FEM and the FEM: 

Ci - + - — +C2— — , iiu(^H\n), 
\v r> \av/ / » V crn/ 



\\u~Uh\\Hi(n) 



< 



Ci{kh)P + C2k{kh)^P, ifp = 0(1) and ||u||^p+i(i^) < kP. 



Here cr is a constant independent of k,h,p, and the penalty parameters. Note that the pollution 
term is 0{k^P+^h^P) if p = 0(1) and < kP and is if w G H^in). 

(ii) The CIP-FEM attains a unique solution for any k > 0, h > 0, p > 0, and 7 > 0. 

(iii) Estimates in the L^-norm are also obtained. 

We remark that the numerical tests in [17] for the linear CIP-FEM show that the penalty parameter may 
be tuned to greatly reduce the pollution errors. 

Error analysis and dispersion analysis are two main tools to understand numerical behaviors in short 
wave computations. The later one, which is usually performed on structured meshes, estimates the error 
between the wave number k of the continuous problem and some discrete wave number (denoted by w) of 
the numerical scheme [T1[11[2I1I3I1|M1[351I1S1I1S]. In particular, it is shown for the hp-FEM (cf. [T[[55]) 
that 



0(fc2P+l/i2p) iffc/l<l, 

k (ekh\'^P\ , , , ^ ekh 

- if fc/i > 1 and 2p+ 1 > . 

p\ Ap J J 2 



O 



By contrast, our pre-asymptotic error analysis, which works also for unstructured meshes, gives error 
estimates between the exact solution u and the discrete solution Uh- Clearly, our pollution error bounds 
in iJ^-norm coincide with the phase difference |fc — wj as above. 

Our analysis relies on two novel tricks and the profound stability estimates for the continuous problem 
given in [311130]. The first one is a modified duality argument (or Aubin-Nitsche trick). The traditional 
duality argument is a crucial step in asymptotic error analyses of FEM for scattering problems, which is 
usually used to estimate the L^-error of the finite element solution by its if ^-error (cf. [H [23l |35l [SO] |40l 
|42]). Our key idea is to use some special designed elliptic projections in the duality-argument step so that 
we can bound the L^-error of the discrete solution by using the errors of the elliptic projections of the 
exact solution u and thus obtain the L^-error estimate directly by the modified duality argument. This 
helps us to derive pre-asymptotic error estimates for the CIP-FEM and the FEM under the condition 
(|1.9p . The second trick is used to prove the absolute stability of the CIP-FEM. One of the authors and 
Feng [26l [27] have developed an approach to analyze the stability of the interior penalty discontinuous 
Galerkin methods which makes use of the special test function vu '■= Vu^ -{x — xq) (defined element-wise) 
and a local version of the Rellich identity (for the Laplacian) and mimics the stability analysis for the 
PDE solutions given in [551 HOI 132] • Here xn is a point such that the domain is strictly star-shaped 
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with respect to it. But the function Vh is not in the approximation space of CIP-FEM which is the same 
as that of the FEM, i.e., the space of continuous piecewise (mapped) polynomials, and thus, it can not 
server as a test function for the CIP-FEM. We circumvent this difficulty by using the projection of 
Vh onto the approximation space of CIP-FEM instead. By using the so-called "Oswald interpolation" 
[T51 1551 [551 WT\ we may show that the difference between and its projection is estimated by the 
jump of Vh at mesh interfaces, and hence, controlled by using the jump term (|1.8I) in the CIP-FEM. We 
remark that the technique for deriving the stability of the linear CIP-FEM developed in the first part of 
the series [47] does not work for higher order methods because it relies on the fact of Auh — on each 
element. 

The remainder of this paper is organized as follows. The CIP-FEM is introduced in Section [2j Some 
preliminary results, including the stability of the continuous solution, the approximation properties of 
the /ip-finite element space, and estimates of some elliptic projections, are cited or proved in Section [3l 
In Section |4l by using an improved duality argument, some pre-asymptotic error estimates in H^- and 

L^-norms are given for both CIP-FEM and FEM under the condition that (|) "^^ ^ is small enough. In 
Section [5J the CIP-FEM is shown to be stable for any fc > 0, /i > 0, p > 0, and 7 > 0. In Section [Bl 
pre-asymptotic error estimates of the CIP-FEM are proved for A: > 0, /i > 0, p > 0, and 7 > by 
utilizing the error estimates for the elliptic projection, the stability results for the CIP-FEM, and the 
triangle inequality. The proofs of two preliminary lemmas on /ip-approximation properties are given in 
Appendix El 

Throughout the paper, C is used to denote a generic positive constant which is independent of h, p, 
k, /, g, and the penalty parameters. We also use the shorthand notation A < B and B > A ioi the 
inequality A < CB and B > CA. A ~ S is a shorthand notation for the statement A < B and B < A. 
We assume that k >1 since we are considering high-frequency problems. For the ease of presentation, we 
assume that k is constant on the domain fi. We also assume that is a strictly star-shaped domain with 
an analytic boundary. Here "strictly star-shaped" means that there exist a point xn G and a positive 
constant cji depending only on such that 

(1.10) {x - xq) ■ n > cn, Vx £ dn. 

Sure the theory of this paper can be extended to the case of polyhedral domains. The extensions to other 
types of boundary conditions, such as PML absorbing boundary condition (cf. [TTl [lOl [18]) and DtN 
boundary condition (cf. |39j ) . will be addressed in future works. 

2 Formulation of CIP-FEM 

To formulate our CIP-FEM, we first introduce some notation. The standard space, norm and inner 
product notation are adopted. Their definitions can be found in [121 1191 . In particular, (•, ■)q and (•, 
for E C dQ denote the L-^-inner product on complex- valued L^{Q) and L^(S) spaces, respectively. Denote 
by (•,•) := i-rh and (•,•) := 

For the triangulation of we adopt the setting of [321 HO] • The triangulation Th consists of elements 
which are the image of the reference triangle (in two dimensions) or the reference tetrahedron (in three 
dimensions). We do not allow hanging nodes and assume — as is standard — that the element maps of 
elements sharing an edge or a face induce the same parametrization on that edge or face. Additionally, 
the element maps Fk from the reference element K to K £ Th satisfy the assumption 5.1 in [40]. For 
any triangle/tetrahedron K G Th, we define Hk ■— dia.m{K). Similarly, for each edge/face e of K £ Th, 
define he '■= diam(e). Let h — max^KeTh ^k- We define 

Sl set of all interior edges/faces of Th, 

£h ■— set of all boundary edges/faces of Th on F. 

We also define the jump [v] of v on an interior edge/face e = dK n dK' as 

r , I j v\k — v\k' , if the global label of K is bigger, 

\ v\k' — v\k, if the global label of K' is bigger. 
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For every e — dK D dK' E £[, let rie be the unit outward normal to edge/face e of the element K if the 
global label of K is bigger and of the element K' if the other way around. For every e <E Sj^ , let n,, — n 
the unit outward normal to dfl. 

Now we define the "energy" space V and the sesquilinear form a^(-, •) on F x 1/ as follows: 



(2.1) 

where 

(2.2) 



ah{u,v) := {VUjVv) + J{u,v) \fu,vGV, 



J{u,v) := V i7e 



du 



dUe 



dv 



One 



and 7e, e G £^ are nonnegative numbers to be specified later. 

Remark 2.1. (a) The terms in J{u,v) are so-called penalty terms. The penalty parameter in J(u,v) is 
ije ■ So it is a pure imaginary number with positive imaginary part. It turns out that if it is replaced by 
a complex number with positive imaginary part, the ideas of the paper still apply. Here we set their real 
parts to be zero partly because the terms from real parts do not help much (and do not cause any problem 
either) in our theoretical analysis and partly for the ease of presentation. 

(b) Penalizing the jumps of normal derivatives was used early by Douglas and Dupont 1221 for second 
order PDEs and by Babuska and Zldmal JS^ for fourth order PDEs in the context of C'^ finite element 
methods, by Baker [9j for fourth order PDEs and by Arnold for second order parabolic PDEs in the 
context of IPDG methods. 

(c) In this paper we consider the scattering problem with time dependence e"^*, that is, the signs before 
i 's in the Sommerfeld radiation condition (|1.4[) and its first order approximation (|1.2[) are positive. If we 
consider the scattering problem with time dependence e~'"*, that is, the signs before i 's in (II. 4p and 

are negative, then the penalty parameters should be complex numbers with negative imaginary parts. 



It is clear that J{u, v) 
(in|)-([Lll), then 



if u G H'^{fl) and v e V. Therefore, if m e _ff^(fi) is the solution of 



(2.3) ah{u,v)-k^{u,v) + ik{u,v) = {f,v) + {g,v), Vw G K 

Let Vh be the hp-CIP approximation space, that is, 

Vh [vh e H\n) : v^ \k oFk £ Vp{K),yK e %} , 

where Vp{K) denote the set of all polynomials whose degrees do not exceed p on K. Then our CIP-FEMs 
are defined as follows: Find Uh S Vh such that 



(2.4) 



ah{uh,Vh) - k^{uh,Vh) + ik{uh,Vh) = (/, fh) + {g,Vh) , Vw^ G Vh. 



We remark that if the parameters 7e = 0, then the above CIP-FEM becomes the standard FEM. 
The following (semi-)norms on the space V are useful for the subsequent analysis: 



(2.5) 



(2.6) 



\l,h 



\Vv\ 



Eke 



dv 



dUe 



LHe) 



1/2 



Mlh + k\\v\\l2^r)[ 



In the next sections, we shall consider the pre-asymptotic stability and error analysis for the above 
CIP-FEM and the FEM. For the ease of presentation, we assume that 7e ~ 7 for some positive constant 
7 and that hx — h. 
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3 Preliminary lemmas 

In this section, we first recall the stability estimates of the continuous problem from Melenk and Sauter 
[39l |40] . Then we consider the /ip-approximation estimates of the discrete space Vh- These estimates 
improve slightly the estimates in |40] when u € H'^{fl). Finally, we analyze some elliptic projections 
which are crucial for our pre-asymptotic analysis. 



3.1 Stability estimates of the continuous problem 

Let V" stand for derivatives of order n; more precisely, for a function m: 17 — )■ R'^, 17 C Mf^, |V"'u(x)p = 
EaGN„^:|al=« The following lemma (cf. 001 Theorem 4.10]) says that the solution u to the 

continuous problem (jl.ip - (jl.2p can be decomposed into the sum of an elliptic part and an analytic part 
u — U£ + uj( where us is usually non-smooth but the iJ^-bound of U£ is independent of k and uj( is 
oscillatory but the iJ-' -bound of uj^ is available for any integer j > 0. 

Lemma 3.1. The solution u to the problem (|l.ip - (jl.2p can he written as u = ug + Uj^, and satisfies 

(3.1) ||wf||//2(n) -I- k\u£\min) + k'^\\u£\\L^n) < CCf^g, 

(3.2) \uA\min) + k\\uA\\L^{n) < CCf^g, 

(3.3) Vp e No, W'^p+^uaWlhu) < CXPk-' max{p, k}P+^Cf,g. 

Here, A > 1 independent of k and Cf,g :— ||/||l2(si) + ||.9||_ffi/2(r) • 
Remark 3.1. An direct consequence of the above lemma is that 

This estimate was proved in i20\ \3S\ \3EI- 



3.2 Approximation properties 

In this subsection we consider to approximate the solution u to the problem ()l.ip - (|1.2p by finite element 
functions in Vh- We give two types of approximation estimates in both L^-norm and the norm |||'|| 
defined in (|2.6p . The first type of estimates can be applied to smooth solution and gives higher order 
convergences both in h and p. The second type of estimates assumes only regularity and gives first 
order convergences (even for higher order elements), but it uses the decomposition given in Lemma l3. II 
and is more subtle than the first one. The proofs are a little long and will be given in the appendix. 

Lemma 3.2. Let s > 2 and fi = min{p+ 1, s}. Suppose u G H'^{Q). Then there exists Uh G Vh such that 
(3-4) II"- "/i|Il2(o) < — hll//»(n), 

(3-5) |||u - uh\\\ < Ccrr-—^\\'^\\H-(n), 

where C,,,:= + 

The following lemma gives approximation estimates of the solution to the problem (ll.ip - (jl.2p . 

Lemma 3.3. Let u be the solution to the problem (|l.ip - (ll.2p . Then there exist Uh G Vh and a constant 
(T > independent of k, p, h, and the penalty parameters, such that 

(3.6) \\u-uh\\mn) < + 

(3.7) lll"-^.|||<Ce„.(- + i(-)')C/,„ 

\p p\crp// 



where Cf^g := \\f\\L^(n) + ||5llffV2(r) and Corr -.^ [l + j + f 



kh\2 
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Remark 3.2. The above estimates are based on the results from J39l\40^ and a duality argument. Similar 
estimates have been given in J391 \40\l as follows 

Our estimates in this lemma improve a little the above estimate in terms of p. See also Remark \A.l\ 

3.3 Elliptic projections 

In this subsection, we introduce two kind of elliptic projections and estimate the approximation errors 
of them. For any u G V, we define its elliptic projections G Vh and G Vh by the following two 
formulations, respectively. 



(3.8) 



ah{ul,Vh) + ik{u+,v,i) = ay,{u,Vh) +ik{u,Vh) Vw/j € Vh, 



(3.9) ah{vh,Uf^) + ik{vh,Uf^) = ah{vh,u) + ik{vh,u) \/vhGVh. 

Let us take a look at the projections from other points of view. Define 



'-{u,v) := (Vm, Vw) ± X! ^' 



du 



dv 



dnp_ 



Then = a/i, (w, w) = ah{v,u) (cf. (j2.ip ). and clearly, e Vh satisfies 
(3.10) aj^{u^,Vh) iikiu'j^^Vh) ^ a'j^iu^Vh) ±ik{u,Vh) Vw/i G V^.. 

In other words, is an CIP finite element approximation to the solution u of the following (complex- 
valued) Poisson problem; 



du 
dn 



-Au = F in n, 
± iku = "0 on r. 



for some given function F and tp which are determined by u. 

Next, we estimate the errors of u^. From (|3.10p we have the following Galerkin orthogonality: 



(3.11) 



'-{u - u^,Vh) ± ifc(u - uf^^Vh) = yvh e Vh- 



We state the following continuity and coercivity properties for the sesquilinear form a^(-, •). Since they 
follow easily from (|2.1I) - (I2.5|) . so we omit their proofs to save space. 



Lemma 3.4. For any v,w gV , 
(3.12) 



\a^{v,w)\,\a^{w,v)\ < \\v\\i^h\\w\\i^h, 



(3.13) Rea^{v,v)±lma^{v,v) = \\v\\l,^. 

Lemma 3.5. Suppose u is any function in H^{i}). Then there hold the following estimates: 
(3.14) 



\u ~ III < inf |||u - z,i||| , 



(3.15) 

where Car, := (l + 7+f 



zheVh 

|u - u^||L2(f2) < Ccrr- inf |||m-z/i||| 

p Zh&Vh 
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Proof. Wc only prove the estimates for since the proof for follows almost the same procedure. For 
any Zh G Vh, let rj =u — , rih=uj^ — Zh- From rjh + rj = u — Zh and p.lip . we have 

(3.16) a/i(r/, rj) + ifc(r/, 77) = 0^(77, u - z^) + ifc(?7, Zh). 

Applying Lemma 13.41 and p.l6p we can obtain that 

hllvi = R-e 0,1(77,77) +Im ah{ri,T]) 

= Re(a/i(77, 77) +1^(77,77)) + Im(ah(77, 77) +1^(77,77)) - fc(77, 77) 
= Re(a/j(77, u~ Zh) + ik{r], u - Zh)) 

+ Im(a,j(77,M - Zh) + ik{r],u- Zh)) - fc||77||i2(r) 

< C (\\7]\\i,h\\u - Zh\\i,h + A:||77||L2(r)||u - Zh\\L-2{r)) - kllvh^ir)- 

Therefore, 

(3.17) |17?|1L + kMh^r) ^ \W - ^kWlh + k\\u - ^hlli.(r). 
That is, (|334l) holds. 

To show (I3.15p . we use the Nitsche's duality argument (cf. [El US])- Consider the following auxiliary 
problem: 



(3.18) -Aw = 77 in fl 

dw 
dn 



Dw 

(3.19) ikw^O on F. 



It can be shown that w satisfies 

(3.20) Mnnn) < hh^n)- 

As a matter of fact, the PDE theory shows that |w|/j2(f2) < C{n)(^ || Au;||^2(52) + ll''^llffi(si) ) i*^^- ^^'^ 
testing (I3.18P by the conjugate of w and taking the real part and imaginary part imply that < 

Let Wfi e Vfi be defined in Lemma [X^ fwith u replaced by w). From p. lip . 

ah{ri, Wh) + ifc(77, Wh) = 0. 
Testing the conjugated (I3.18[) by rj and using the above orthogonality we get 

\\v\\h{n) ^-{u~ ■«+, Aw) = ah{u - m+, w) + ik{u - u+, w) 
= 0,1(77, w - Wh) + ifc(r/, w - Wh) 

< ||?7||i,;i||w - Wh\\i,h + fc||7/||i2(r)||u; - WhWh^iv) 

< III77III \\\w - Wh\\\ 

^ llhlll Ccn--\\w\\H2(n), 

p 

which together with p.l4p and p.20p gives p.lSp . The proof is completed. □ 



By combining Lemma [33] and Lemma l3.3l we have the following lemma which gives the error estimates 
between the solution to the problem (|l.ip - ()1.2p and its elliptic projections. 

Lemma 3.6. Let u be the solution to the problem (|l.ip - ()1.2p . Then there hold the following estimates: 



/i^ h (kh\v\ 



(3-22) ll"--"ll-^(-)^^-^-(^ + ^(S)>/- 



where Ccrr := (1 + 7 + 



kh\2 
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4 Pre- asymptotic error estimates by duality argument 

One crucial step in asymptotic error analyses of FEM for scattering problems is performing the duality 
argument (or Aubin-Nitsche trick) (cf. [H [23l [35l [39l |40l [42] ) . This argument is usually used to estimate 
the L^-error of the finite element solution by its ff^-error. In this section we use a modified duality 

argument to derive pre-asymptotic error estimates under the condition that (|)''^^-^ is sufficiently 
small. The key idea is to used the elliptic projections from the previous section in the duality-argument 
step so that we can bound the L^-error of the discrete solution by using the errors of the elliptic projections 
of the exact solution u and thus obtain the i^-error estimate directly by the modified duality argument. 

Theorem 4.1. Letu anduh denote the solutions of (|l.ip - (|I.2p and , respectively. Suppose < 7 ^ 1- 
Then there exist constants Cq and a > independent of k,h,p, and the penalty parameters, such that if 

then the following error estimates hold: 

(4.2) < fl + inf \\\u-zh\\\, 

V p \apJ J zueVh 

(4.3) \W -Uh\\L^{n) <(- + -(— Y) inf \\\u - Zh\\\ . 

\p p \ apJ I zh£Vh 

Proof. First, we estimate the L^-error by introducing the dual problem and using the elliptic projections of 
the solution to the original continuous problem and of the solution to the dual problem. Let Ch := u — Uh- 
Consider the following dual problem: 

(4.4) -Aw - fc^w = Ch in fi, 

(4.5) - ifcw on F. 

an 

Let uj^ be the elliptic projection of u defined by (|3.8p and let G Vu be the elliptic projection defined 
as (13. 9p . From (|2.3p and (12.41) we have the following orthogonality, 

(4.6) ah(ch,Vh) - k^{eh,Vh) + ik{eh,Vh) ^Q, Vw/i G V/j. 
Testing the conjugated (14.41) by e/,. and using (|4.6p with Vh — w^, (|3.8p - p.9p we get 

(4.7) l|e/Jli2(o) = (Ve,,, Vw) - k^{ch,w) + ik{eh, w) 

= Ohich, w) - k^{eh,w) + ik{eh, w) 
= ah{eh,w~ wy^) + ik{eh,w- w^) - k'^{eh,'w- w^) 
= ah{u - uj^,w - w'f^) + ik{u - u'l,w - w^) - k'^{eh,w - w,7). 
Similar to Lemma 13.61 we may show that 

1 



(4.8) \h-Wy\\\<Corr(- + -( 

III III \p p\ 



p\ap 
2 //i^ h /kh\P 



Ch 



(4-9) \\^-<\\L^n)<C^r.[-^ + ^[^)' jWekh^^a) 



From (|47l) . Lemma [331 Lemma [331 and ([4 ^ - ([49|) . we have 

WehWh^n) < h - u+\\i^h\\w - w-\\i^h + k\\u - u+\\^^^^^ \\w - w' \\ 

\^h\\L^(n) \\w - ||^2(n) 



k^ 



< CC^J- + -(—)") \\eh\\L2u-^. inf 



\U — Zh 
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Noting that Corr = 1 + 7+ ^ ^ l + there exists a constant Cq independent of k,h, p, and the penalty 
parameters such that 



and as a consequence, 

^ fh 1 fkh\P\ 

(4.11) ll^''ll-^(-)^S(p+pW )..¥.J1I"-^'^II1- 



That is, (1131) holds. 

Next we turn to prove (14. 2p . Let e V/i be the eUiptic projections of u defined by p.8p - (l3.9p and 
denote by = w — . It follows from Lemma 13.41 and (|4.6|) that 

\\(^h\\\,h =R-eah(e;i,eft) + Ima,i(eft, e^) 

= Re{ah{eh, en) - k'^{eh, eh) + ik{eh, eh)) + k'^{eh, eh) 

+ Im(aft(e,j,e;i) - k^(eh,eh) + ik{eh,eh)) - fc||eft|||2(r) 
= Re{ahieh,C) ~ ^^(eh, C) + ifc(eh, O) + fc^||e^|||2(o) 

+ Im(aft(eh,C") - fc^(eh,C") + ifc(e/i,C")) - fc||eh||i2(r) 
= Re(a;,(C+,r) - fc'(e,„r) +ifc(C+,r» + fc'||e„||i.(n) 

+ Im(a^(C+, D - k^{eh, C) + ifc(C+, C » - k\\eh\\h^r) 
<2||C+l|i,/.||r 111,/. + fc'lirili2(n) + 2fc||C+llL^(r)lir llL^(r) + 2k^eh\\h^n) " k\H\\h(ry 

Therefore, from Lemma [5751 and noting that ^,7 J; 1, 

ll|e/i||| < inf |||u - z/illl + A:||eh||L2(Q). 

which together with (14.31) implies (14.21) . This completes the proof of the theorem. □ 

Remark 4.1. (i) Noting that j/7 = then the CIP-FEM becomes the FEM, the above theorem and the 
three corollaries below hold for the standard FEM. 

(ii) Noting that the wave length is the condition (|4.ip roughly says that about ^{^)''~'^^ degrees 
of freedom are needed in one wave length. 

(iii) From (|4.10p , clearly, the theorem holds under the following condition which is a little more general 
(and a little more complicated as well) than (14. ip ; 

, , fc^/i^ k (kh\P+^ _ 

(4.12) — - — I — I — IS sufficiently small. 

p^ p\apJ 

(iv) // 

kh k (kh\P 

(4.13) \ — I — is sufficiently small, 

p p\apJ 

then CIP-FEM, as well as the FEM, is pollution free in the norm \\\-\\\, i.e., \\\u — Uh\\\ ^ infz^gy,^ |||u — z/i||| . 
The condition (|4.13p improves a little the condition of ^ + fc(^)'' small enough given in I39j for the 

FEM, because, for example, if we choose k{^y — c for some fixed constant c so that the latter condition 



does not hold, then ^ + ^{^)^ — (f ) + ^ may be small enough for an appropriate choice of p and 
large k. 

(v) The pre- asymptotic error estimates for the linear CIP-FEM and FEM (p — 1) in two and three 
dimensions are given in the first part of this series f^Tf . For the pre- asymptotic error estimates for the 
linear and hp version of the FEM in one dimension, we refer to \34^ and J35f . respectively. To the best 
of the authors' knowledge, there have been no pre- asymptotic error estimates for the higher order FEM 
and CIP-FEM (p> 2) in two and three dimensions given in the literature so far. 
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From Remark I4.1f iv) we have the foUowing corollary which gives practical sufficient conditions for 
asymptotic estimates of the CIP-FEM and the FEM. 

Corollary 4.1. Let u and denote the solutions of (|l.l|) - (jl.2p and ^2.4^ .respectively. Suppose < 7 < 
1 . Then the following asymptotic error estimate 

III" - ■"'illl ^ inf |||u - Zh\\\ 
zheVh 

holds under any of the following three conditions: 

(4.14) ( — J — is sufficiently small, 

\pJ p 

(4.15) p > (lnfc)2 and {— \ — is sufficiently small, 

\pJ p 

k h 

(4.16) p> Ink and is sufficiently small. 

P 

Proof. (|4.14p impUes (I4.13|) and the conditions (|4.15p - (|4.16l) are all sufficient conditions of (|4.14l) . □ 

Remark 4.2. (i) Our asymptotic estimates hold for both the CIP-FEM and the FEM. Melenk and Sauter 
\39\. \40^ proved recently the asymptotic error estimates for the FEM under either of the condition (I4.16|) 
and the following one: 

(4.17) p = 0(1) fixed independent of k and k^~^ph is sufficiently small. 

The above condition is similar to the condition (|4.14p but we prefer the later one because it does not 
require p = 0{\). 

(ii) The condition (j4.15l) says that the theorem 14.11 gives actually asymptotic estimates in the norm 
lll'lll under the further condition of p > (lnA:)2 . This condition is not strict from the practical point of 

view. For example if we take p ~ {\nk)^ then p grows slowly as k increases, and so is (-j''^^ — 
which is o(k^) for any positive constant e. 

(iii) The condition ()4.14|) is convenient for fixed p . The condition (|4.16|) requires p > In/c and is perfect 
when implementations of sufficiently high order CIP-FEM or FEM are available. And the condition (|4.15|) 
which requires merely p > (lnfc)2 is suitable when only "medium" order methods can be used. 

From Theorem 14.11 and Lemma 13. 3[ we have the following corollary which gives estimates for 
regular solutions. 

Corollary 4.2. Suppose the solution u G II^{^). Under the conditions of Theorem there hold the 
following estimates: 



(4-18) <[- + -[-] +-J - ]C. 



'ti I /kh\P k /kh\'^P 
(4.19) u,,\\i^.^a) < i'^ + ^(—YlCf, 



p p\ap/ p^ \ap/ 

1 /fc/l\2p 

— + — — 

p^ V ap 

Remark 4.3. (i) Ifp = o((lnfc)2 j then the condition (j4.ip is weaker then the condition (j4.14p for asymp- 
totic error estimates, and the estimate (|4.18p is a pre- asymptotic one with pollution term O^p-(^)^^^ . 

(ii) Pre- asymptotic error analysis and dispersion analysis are two main tools to understand numerical 
behaviors in short wave computations. The later one which is usually performed on structured meshes 
estimates the error between the wave number k of the continuous problem and some discrete wave number 
m\M\M\M\3^\M\4^\45^- In particular, it is shown for the hp-FEM (cf. \3Bi) that 

0{k^P+^h^P) ifkh-^l, 
^/ k / ekh\'^P\ , ekh 

) ^fkh»land2p+l>—. 
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By contrast, our pre- asymptotic error analysis gives the error between the exact solution u and the discrete 
solution Uh and works for unstructured meshes. Clearly, our pollution error bounds in H^-norm coincide 
with the phase difference \k — lj\ as above. 

(iii) Our pre- asymptotic estimates in Theorem \4. 1\ and Corollary \4.2\\J73\ hold in particular under the 
following condition 

(4.20) p — 0{1) fixed independent of k and k^'^p+^^h is small enough, 

which extends the asymptotic range given by (|4.17p . 

By combining Theorem 14.11 and Lemma 13.21 we have the following corollary which gives estimates for 
H'^ regular solutions (s > 2). 

Corollary 4.3. Suppose the solution u e H'^{Q), s > 2. Let ^ — min{p+ 1, s}. Under the conditions of 
Theorem \4.1\ there hold the following estimates: 



(4.21) <(! + _(_) 

II I, ^ fh 1 / kh\P\ h'^^'^ , 

(4.22) <[- + -[-) )^\\u\\HHn), 



where the invisible constants are independent on p, h, k, and the penalty parameters, but may depend on 
s. In particular, if u E H'p^^(VL) is an oscillating solution in the sense of fSSI Definition 3.2], i.e. 



l"lli/p+i(o) ^ 



then the following estimates hold: 

^ „, .nkh\P k/kh\'^P\ 
(4.23) \\-u.h.^Cip){{-) +-(-) ), 

I, I, ^ C(p) / i'kh\P+^ kfkh\2p\ 



Remark 4.4. (i) Corollary \4.3\ certainly holds for s = 2. We exclude the the case s — 2 in this 
corollary because the estimates in (I4.2ip - (|4.22p with s = 2 are worse than those in (|4.18p - (|4.19p . For 
s > 2, Corollary \4.S\ gives optimal convergence order in h and p while Corollary \4.2\ gives only first order 
convergence in h and p. 

(ii) If p = 0{1) and u is an oscillating solution then (j4.23p shows clearly that the pollution error in 
H^-norm is again bounded by O^k^^^^h^^). 



By combining Lemma lOI and Corollarv l4.2l we have the following stability estimates for the CIP-FEM 
(FEM). 

Corollary 4.4. Suppose the solution u G H^(fl). Under the conditions of Theorem \4.1\ there hold the 
following estimates: 

IWhWij, + fc|lw/illL2(o) < 

and hence the CIP-FEM is well-posed. 

Remark 4.5. If^~0, then this corollary gives the following stability estimate for the standard FEM 
under the same condition. 

\\^Uh\\L2^n} + k Wuhh^n) < Q.s- 
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5 Stability estimates for the CIP-FEM 

In the previous section we have shown that the CIP-FEM and the FEM are stable under the condition 
that ^ < Co(f ) ''^^ and < 7 < 1. The goal of this section is to derive stability estimates (or a priori 
estimates) for the CIP-FEM ((2^ for any k,h,p, and 7 > 0. 

Some stability estimates have been proved for the discontinuous Galerkin methods (cf. [26l [27l [28] ) 
and the spectral-Galerkin methods (cf. [44 ). Their analyses mimic the stability analysis for the PDE 
solutions given in [301 1311 [3H]- The key idea is to use the test functions v = Uh and v = (x — xq) ■ Vut, 
respectively, and use the Rellich identity (for the Laplacian), where xn is a point such that the domain 
n is strictly star-shaped with respect to it (see (jl.lOp ). As for our CIP-FEM (|2.4p . although the test 
function Vh = Uh can be still be used, the test function Vh = {x — xq) ■ Vu^i does not apply since it is 
discontinuous and hence not in the test space Vh- We surround this difficulty by taking the projection 
oi {x ~ xn) ■ '^Uh to the finite element space Vh- To analyze the error of the projection, we recall the 
so-called Oswald interpolation operator Jqs, which has been analyzed in [Tni|331IMlllI] • Let Wh be the 
space of (discontinuous) piecewise mapped polynomials: 

(5.1) Wh := n KiFK'iK))- 

Lemma 5.1. There exist an operator Iqs ■ Wh 1— > Vh and a constant Cp depending only on p, such that, 
for all K G Th, the following estimate holds: 

(5.2) yVheWh, Wvh-^OsVhWL-^iK) <Cphl^ ^ IIK]||L2(e), 

where Cp — p^^ for Cartesian meshes and Cp = p^~ for triangulations {d = 2, 3). 

Lemma 5.2. Let Qh be the projection onto Vh and let Cp he the constant in Lemma \5-1[ Then the 
following estimate holds: 

(5.3) yvhEWh, lk/.-0/.«/.||L^(r.)<C^p(E'^-IIMIIi^(e))'' 



2 / 

(5.4) wvh e Wh, \\vivh - QhVh)\\mn) ;$ ^ E '^-IIK: 



H 
h 



2 

\L^(e) 



Proof- (15.31) follows from \\vh — QhVh\\L'^{Th) — ~ ^OsVh\\L^{Th) ^^^'^ Lemma [5.11 (|5.4p follows from 
the inverse inequality and (|5.3D . □ 



We cite the following lemma [26l Lemma 4.1], which establishes two integral identities and play a 
crucial role in our analysis. 

Lemma 5.3. Let a{x) := x — xq , v £ Yi H'^i^)' K €Th and eg £^_^. Then there hold 

Ken 

(5.5) d\\v\\l2,K) + '2^<v,a-Vv)K ^ / a • n^lw p, 

JdK 

(5.6) {d-2)\\yv\\l2(K)+2Re{\/v,\7{a-\7v))K ^ f a-niflVvj^. 

JdK 

Here xq is a point such that the domain is strictly star-shaped with respect to it (see (jl.lOp ]. 
Remark 5.1. The identity (j5.6p can be viewed as a local version of the Rellich identity for the Laplacian 

A (cf imj)- 
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Lemma 5.4. Let E Vh solve (|2.4p . Then 



(5.7) \\^Uh\\L2(n) - k^hhWh^n) < 2 11/11^2(1^) ||w/J|i2(f^) -r - iiyiiL2(r) 



(5.8) 



Ehe 



duh 



One 



L2(e) 



fc|kh|li2(r) < 2 11/11^2(0) ll"/»llL2(n) + Il5lli2(r) 



Proof. Taking — in (I2.4p yields 

(5.9) ah{uh, Uh) - k^\\uh\\'i2(^) + ik\\uh\\\2(Y) = (/, "ft) + (ff, "/i) 

Therefore, taking real part and imaginary part of the above equation, we get 
(5.10) 

(5.11) 



2 



Ehe 



9ne 



From (f5ll|) . 



Ehe 

P 



duh 



due 



1 fc 



which implies (|5.8p . From (j5.10p and 



1 fc 

1 2 

< 2||/|lL2(i^) lk/»llL2(o) + T ll5llL=(r) • 



That is, (15. 7p holds. This completes the proof of the lemma. 



□ 



From ((5?7)) and ((O)) we can bound || Vuft||^2(f2) ^^'^ jumps of across e G In order to get 
the desired a priori estimates, we need to derive a reverse inequality whose coefficients can be controlled. 
Such a reverse inequality, which is often difficult to get under practical mesh constraints, and stability 
estimates for Uh will be derived next. 

Theorem 5.1. Let ut G Vt solve (|2.4p and suppose 7e ~ 7, hK,he ~ h. Then 
(5.12) k\\uh\\mn) + \\uh\\i,h < CstaM(/,3), 

w/iere 



(5.13) 
(5.14) 



M{f,g) |l/|lL2(^,) + ||g|U2(r), 

7/12 



Proof. We divide the proof into three steps. 

Step 1: Derivation of a representation identity for \\uh\\L^(n)- Define Vh by Vh\K = Oi ■ 'S/uh\K for 
every K ETh, denote by Wh = QhVh, hence Wh € V/j. Using this Wh as a test function in (|2.4p and taking 
the real part of resulted equation we get 



(5.15) 



k'^ Re{uh, Wh) = Re((/, Wh) + (g, Wh) ~ ah{uh,Wh) ~ w/i))- 
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It follows from ((53|) . (ICTI) . and ((5?T5|) that 

(5.16) 2e\\un\\l.^^)^k^ I c.-nK\un\'^{d^2)e\\un\\l.^^^ 



~ 2fc^ Re(u/i, Wh) - 2fc^ Re(u/i, w/j - Wh) 
= fc^ ^ / a • + (d- 2)Re((/,M;i) + - a,i(u,i,u,i)) 

+ 2 Re((/, Wh) + (5, w^/i) - ah[uh, Wh) - ik{uh, Wh)) 
= k^ [ a-nK\uh\^ + {d-2)Re{{f,Uh) + {g,Uh)) 

+ 2Ke{{f,Wh) + {g,Wh))- (('^-2)||Vu,,|| 

KeTh 

- 2ReJ{uh,Wh) +2klm{uh,Wh). 
Since Uh is continuous, we have 

(5.17) Y / a • = 2 ^ Re(a • ne{Mft}, + (a • nn, 

= {a ■ nn, \uh\'^)- 

Using the identity |ap — |6p = Re(a + fe)(a — b) followed by the Rellich identity (|5.6p we get 

(5.18) ^ ((d-2)||VMft|| 

h{K) + 2 Re(Vu/i, Vw/i)_ff ) 
= 5] ((d-2)||V7.,.|| 

L2(x) + 2 Re(Vu/i, VvhjK) 
+ 2Y Re(Vw,„V(u;ft- 



^Y a-nK\Vuh\^ + 2Y ^'^^'^^h,'^{wh-Vh))K 
KeTh •'^^ KeTh 

= 2 ^ Re(a • n^{\Juh}, [Vw^De + (a • no, iVu^p) 



+ 2 ^ Re(Vw,„V(w,, - w,0)k. 
KeTh 

Plugging ([CTT)) and (|5TTS|) into (ISH)) gives 

(5.19) 2fc2||M^||i.(,y = (d- 2)Re((/,M^) + (g, u;,)) + 2 Re((/, i>;,) + {g,Vh)) 

+ fc^(a • no, \uh\^) + 2fclm(uft,w,i) - (a • no, \Vuh\'^) 

- 2 ^ Re(a • ne{Vu/i}, \^Uh])e - 2ReJ{uh,Vh) 
^^£'h 

-2ReJ{uh,Wh~Vh)~2 Y 'Re{'Vuh,'V{wh - Vh))K 

KeTh 

+ 2klm{uh,Wh - Vh) + 2Re((/, w/i - Vh) + {g,Wh - Vh))- 

Step 2: Derivation of a reverse inequality. We bound each term on the right-hand side of (|5.19p . We 
have 

(5.20) 2Re((/,«,0 + {g,Vh)) <C\\f\\h^,r) + I II V^,.||i.(o) 

+ C\\g\\Ur) + "fY l|V^''lli^(e)- 
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It is clear that 

(5.21) k^{a ■ no, \uh\^) < CeWuhWh^^r)- 
It follows from the star-shaped assumption on that 

(5.22) 2fclm(u,„t;,,) - (a-n^jVyftp) 



< Cfc ^ ||M,i||i2(e)||Vu;i||i2(e) - Co ^ ||Vu,i|| 



2 



2 

L2(e 



For an edge/face e G £[, let and denote the two elements in Th that share e. We have 



1 



(5.23) 



-2 E Re{a ■n,{Wuh},[^Uh])e < C J2 {j^)' W^^hh^K^U 

<^l|V«.|li.(o)+C^i^^E 



duh 



dn, 
duh 



L2(e 
2 



From ()2.2p . the trace and inverse inequalities (cf. Lemma lA.ip . we have 



(5.24) 



-2ReJ{uh,Vh) = -2Re i^e^ 



duf, 



duh 



< 



E 



ee£! 



duh 



< 



Ehe 



duh 



dUe 

4 



dvh 



dn, 



,2. i 



K=K^,K' 



n2. 1 



K=Kc.K' 



^ 1 IIY7 l|2 , TP 



he 
pZ 



duh 
One 



Next we estimate the terms containing Vh — Wh — Vh ~ QhVh- From Lemma 15.21 and the definition of Vh-, 



(5.25) 
(5.26) 



duh 



72/1 



dne 

duh 
dne 



L2(e) 
2 
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From (|5.26p and the trace inequality in Lemma [A. 11 



(5.27) 



-2ReJ{uh,Wh-Vh) < 



duh 

dUe 

duh 
dur. 



||[V(wh - Vh)]\\L^e) 



(-^)'l|V(wh-'y/i)||L2(K) 



duh 



< C 



Ehe 



dn. 

duh 

due 



L2(e 
2 



\Vh~Wh\HHn) 



And we have, 



(5.28) -2 ^ Re(Vu,„V(w,,-i;,,))K < 2 ^ ||Vu,,||i2(^)||V(w,. - w,0|li2(x) 

A'eTi, A'eTft 

1 1,2 , ^ I ,2 



< 



1 



L2(r2) 



7/l2 



dUe 



From (|5.25p and Lemma [KH 

(5.29) 2fclm(w;i, W,^ - W/j) < Cfc ^ ll'"'i|lL2(e) II"''' ~ "'«llL2(e) 



2v i 



< Ck ||M/i||i2(r) (^) ' llwft - Vh\\L2^j-^) 



^»iKii™(?)*f (i: 



he 



k \\uh 



1 2 \ ^ 

li^(r) + Te-^ 



L2(e) 
2 



L2(e) 



From ([05)) . ([QS]) . and Lemma [Ol 

(5.30) 2Re((/,w;,,-Wft) + (g,w;ft-i;,,)) 



< 



< 



ll/llL2(a) ||u;/i - vh\\L2^r^) + ||5llL2(r) ( IT ) ' II""* " ""^h^n) 



( ll/llL2(n) + ll5llL2(r) f J2 >!5 



<II/Ili2(n) + 11.91'' 



(7ft) 



9ne 



L2(e) 



L2(r) 



E 



duh 



due 



Plugging (|Og)) ~ (IOi)) and (j^TTfl - ljOgi) into (jgl^ we obtain 

2fc2||u,||i.(^) < K/,^.,,) + {g,u,,)\ + C(||/|li.(^) + ||g||i.(r)) 

+ I l|Vuh|li2(n) - ^ 



E l|Vu/.|li2(,)+c(fc+p2Cp( 



c 



^7ft 



k\\uh"^ 



7e — 



duh 

dUe 



L^^) 



L2(e) 
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From (|5.8p and the facts 



we have 



Step 3: Finishing up. It follows from (|5.7p . (|5.13p . (|5.14p . and the above inequality that 
< - 2 + C'(II/IIl2(s2) + Il5lli2(r)) 

+ CfcCsta( ll/llL2(f2) ||"/i|Il2(o) + T ll3llL2(r) ) 



3fc2„ „o 1 



which together with (|5.8p implies (|5.12p . The proof is completed. □ 



Since scheme (j2.4l) is a linear complex-valued system, an immediate consequence of the stability 
estimates is the following well-posedness theorem for (|2.4I) . 

Theorem 5.2. T/ie CIP-FEM (|2.4I) /las a unique solution for k > 0, h > 0, p > 1 and 7 > 0. 

Remark 5.2. (i) For f/ie general case when the meshes may be nonuniform, Theorem \5.1\ and Theorem \5.2\ 
still hold with h replaced by h = min^gT-^ hx- The proof is similar and is omitted. 

(ii) When p — I, a better stability estimate is derived in |^7| / by another approach taking advantage of 
Auh = on each element. 

(iii) The stability may be enhanced by an over-penalized technique, just like what is done for the hp- 
IPDG method (cf. \27l ). More precisely, if we replace ah{uh,Vh) in the CIP-FEM (|2.4p by the following 
sesquilinear 



1 fhe^'^'~^ 
al{uh,Vh) = (Vu,j,Vw/i) + E E ( '7^ 



p- 



d^Uh 




d^Vh 




1 


. dni _ 



with 7j^e — Ij > O7 J = lj2, • • • ,q, '2 < q < p, then, after some tedious but similar derivations, we may 
show that (j5.12p holds with Csta replaced by the following stability constant: 



k\ ^ih (7172)2/1 h ^^71^ /i i<j<g-i V7J+1 / h 



4g-4 

2(j-l) p'^C^''^^ 

We omit the details here. Clearly, if we choose 7j ~ {Cph) 29-1 ^ 1 < j < q, then C^^^ < ^ — , 

kh^^ 

where Cp is defined in Lemma \5.1\ Recall that the best stability constant obtained so far for the hp-IPDG 



method is O(^) (cf ^ Theroem 3.3]). 
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6 Pre- asymptotic error estimates for the CIP-FEM by using the 
stabihty 

In this subsection we shall derive error estimates for scheme ()2.4p . This will be done by exploiting the 

linearity of the Helmholtz equation and making use of the stability estimates derived in Theorem 15 . 1 1 and 

the projection error estimates established in Lemma [3.61 

Let u and Uh denote the solutions of (|l.ll) - (ll.2p and (|2.4p . respectively. Recall that Ch ^ u — Uh- Let 
be the elliptic projection of u as defined in (|3.8I) . Write Ch — rj — with i] := u — u^, ^ := Uh — uj^. 

From the Galerkin orthogonality ()4.6|) and p.8|) we get 

(6.1) ah(^, Vh) - k^i^, Vh) + ifc(^, w/i)r = 0^(77, Vh) - k'^ir], Vh) + ik{r], Vh)T 

= -k'^{r],Vh) yvheVh- 

The above equation implies that ^ S V^i is the solution of scheme (|2.4p with sources terms / = —k'^rj and 
g = 0. Then an application of Theorem 15.11 and Lemma 13.51 immediatelv gives the following lemma. 

Lemma 6.1. ^ = — satisfies the following estimate: 

(6.2) m\\mn) + m\iM<Cst.C,„— inf 

p z^eVh 

We are ready to state our error estimate results for scheme (|2.4I) , which follows from Lemma 13.51 
Lemma 16.11 and an application of the triangle inequality. 

Theorem 6.1. Let u and Uh denote the solution of (|l.ll) - (|1.2p and (|2.4p . respectively. Suppose u G 
H%n),s > 2. Then 

( k^h\ 

(6.3) ||u - U/illl./i < (1 + CcrrCsta ) inf |1 I'M - Z,i ||| , 

V p ) Zh&Vh 

(6.4) ||w-W/J|L2(m <Corrfl + fcCsta)- mf 

where Con- := (l + 7+ -y-)^ and Csta is defined in Theorem \5.1\ 

Remark 6.1. (i) For the linear CIP-FEM, we have shown that J^7| /, if ^ \, then the following pre- 
asymptotic error estimate hold for any k,h > 0: 

\\u~Uh\\i,h < {Cikh + C2mm{eh\l})Cf^g. 

(ii) The theorem can also be extended to the over-penalized CIP-FEM. In particular (cf. Remark \5.2\ 

(Hi)), if p > q > 2, and jj ~ {C^h) ^-j-i , 1 < J < q, then the following pre- asymptotic error estimate 
holds for any k,h,p > 0; 

/ kh\h t-i-^ -1 
\\u~Uh\\i,h<Ci inf z/illl +C2 1 + — kp^Cp'-^h-^ inf 

Zh<^Vh \ p / Zh^Vh 

The details will be reported in a separate work. 

A Approximation by hp-finite element 

In this appendix, we prove Lemma 13.21 and Lemma 13.31 



A.l Proof of Lemma [3.21 

We need the following trace and inverse inequalities (cf. [121 [IS])- 
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Lemma A.l. For any K eTh and z eVp{Fj^^{K)), 

We recall the following well-known hp approximation properties (of. [71 [501 139) ): 

Lemma A. 2. Let fi — min{p + 1, s}. Suppose u e H'^{Vl). 
• There exists Uh G Wh such that. 



(A.l) \\u-Uh\\H^(rH)- 

• There exists Uh € Vh such that 
(A.2) 



Ken 



h'^~^ 

<^MHHn), J = 0,1, 2. 



ht"-^ 

pS 



u- UhllnHn) ^ -—-\\u\\H-in), j = OA- 



Here the invisible constants in the two inequalities above depend on s but are independent of k, h, and p. 

Clearly, (13.41) holds (cf. (|A.2I) '). It remains to prove p.Sp . It follows from the (|A.2p and the trace 
inequality that 



(A.3) 

From the local trace inequality 



1 1 /l'""2 



(A.4) 
we have 



E E 

KeTh e£dK 



d{u ~ Uh) 



drip. 



< h- 



On the other hand, from Lemma lA.ll 



< \\u - uhWH^TH) + t;"!!"'' ~ "^hWH^n) 



Therefore, from (|A.1I) and (|A.2I) . 



E E 



d{u - Uh) 



due 



< 



p- 



^ p2s-2 W^WH'in)- 



Then it follows from (|X2|) and (|2?5|) that 



\u - Uh\\l,h ^ I|V(m - Uh)\\l2(^a) + E >^ 



g(M - Uh) 

One 



By combining the above estimate, (jA.3| . and (|2.6p we conclude that p.5p holds. This completes the 
proof of Lemma [ 
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A.2 Proof of Lemma ISTSl 

Lemma IA.2I can also be used to approximate the elliptic part of the solution u. The following lemma 
gives an approximation to the analytic part of the solution u. 



< 



Lemma A. 3. Suppose w_4 is analytic and satisfies (I3.2|) - ()3.3p . 

• There exists UA,h G Wh such that 

(A.5) \\uA-UA,h\\H^{TH)-^{^ \\uA-UAJifH2(^K) 

KeTh 

• There exists iiA.h G Vh such that 

(A-6) \\uA~UA.H\\HHn)<{^ + ^{-) \Cf,„ j ^ 0, 1. 

Here a > is some constant independent ofk,h, and p. 

Remark A.l. The following estimate is proved in J39[ Theorem 5.5] (see also 14 0[ Proposition 5.3]): 

, \ ^ ( kh\fh fkh\P 

^i^^^[\\y{uA-Vh)\\L2^n)+k\\uA-Vh\\L^n))l^ V^^)\p + \-^) 

Clearly, a combination of our -estimate and L"^ -estimate in (jA.6|) gives better upper bound than the 
above estimate. This is because our Lp' -estimate improves that implied in the proof of '^39, Theorem 5.5]. 

Proof. Following the proof of [39l Theorem 5.5], we start by defining for each element K £ Th the constant 
Ck by 

(A 7) c^. V IIV^^^IIW) 

^ ■ > J^^(2Amax{p,fc})2p 

and we note 

(A.8) WyuAhHK) < {2Xmax{p,k})PCK Vp € No, 

KeTh 

Then, |391 Lemma C.2] and a scaling argument provides an approximant UA,h G Wh which satisfies for 
J =0,1, 2, ^1 



Summation over all elements K £Th gives 



KeTh 



ap/ 



KeTh 



The combination of the above inequality and (jA.9|) yields 

(A.10) \\uA~UAAHHTh,<[-,h ij^) +-r(J J^/^' 



Pre- asymptotic Analysis of CIP-FEM and FEM 



22 



Furthermore, we estimate using h < diam 17 and a > (independent of h) 



(A.ll) 



we therefore arrive at 



h \P+i hP-^ o , h^-^ 

TP < 



h + a 



p2-i {h + cr)P+l 



< 



/i^"^ / diam V, \p-^ 



Vcr + diam 



p^-j" V ft, + cr 



3-1 r)2 



P 



2-j 



< 



1 ft2-j h^-j ^khy 



Cf, 



-k i p \apJ 
That is, (|X5|) holds. 

Similarly, |39l Lemma C.3] and a scaling argument provides an approximant uj^^h G Vh which satisfies 
for j — 1, 



Summation over all elements K Cz Th gives 



\\uA ~ UA.hWn^n) 



< 



h \2p+2 k"^ ^khyp 



p^ \ap 



Ken 



The combination of the above inequality and (|A.9p , (|A.lip yields 



< 



h 1 /kh\P 
hp p\apJ 



C 



We introduce the sesquilinear form a(-, •) := (Vu, Vw) + (u, w) and define the elliptic projection UA,h G ^^/i 
of UA by 



(A.12) 
Then, 



Therefore, 
(A.13) 



a{uA,h,Vh) ^ a{uA,Vh) Vw^ G 14. 



- ""-4,^11^1(0) = a{uA ~ UA,h, UA - UA.h) = a{uA - UA,h, UA - UA,h) 



< 



\uA - UA,h\\m{n)\\uA ~ UA.h\\H^{n)- 



\uA - UA.hWHHn) < Ih^ - UA,h\\m{n) < 



h 1 /kh\P 
kp p\<jpJ 



C 



f-g- 



-Aw -\- w — UA — u-A.h in 
dw 



dn 



= 



on r. 



To show jXeJ, we use the Nitsche's duality argument and consider the following auxiliary problem: 
(A.14) 
(A.15) 

It can be shown that w satisfies 

(A.16) \\uj\\m(n) ^ \\uA - UA,h\\L2{n)- 

Let Wh G Vh be defined in Lemma [A. 2 1 Testing the conjugated (jA.14[) by ua — UA,h we get 

\\uA - UA,h\\mn) = {uA - UAji, -Aw + w) 

= (V(u^ - UA.h),^Uj) + {uA - UA,h,w) 
= a{uA - UAji, w) = a{uA - UA,h,w - Wh) 



< 
< 



\UA - UA.h\\H^n)\\w - WhllH^n) 
\uA - UA.h\\H^n)-\\w\\H^{n), 
which together with ()A.16[) and (|A.13p gives (|A.6p . The proof is completed. 



□ 
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Next we prove Lemma 13.31 We consider to approximate U£ and u^, respectively. First, from 
Lemma lA. 21 there are two functions iis^h G Vh and U£_h & Wh such that 

(A.17) \\u£~us,h\\Hi{n)^^^Cf.g, j = 0, 1, 

P 

(A.18) W-n£A\HHTu)1^^Cf,g, j = 0,l,2. 

P ■' 

On the other hand, from Lemma [A. 3[ there exist UA,h G Vh and UA,h G Wh satisfying (jA.6p and (lA.sp 
respectively. Let Uh — u^^h + u^^/i, = Uf^/j + uj^^u- It follows from the triangle inequality that 

(A-20) ||,_^,||^,(^^^<(^ + /^0^)C,,, j=0,l,2. 

Therefore, p.6p holds. 

p.7p can be proved by following the same procedure as that for deriving p.5p . We omit the details. 
This completes the proof of Lemma 13.31 

Acknowledgments. The author would like to thank Professor J.M. Melenk for helpful communica- 
tions. 



References 

[1] M AiNSWORTH, Discrete dispersion relation for hp-version finite element approximation at high wave 
number, SIAM J. Numer. Anal., 42 (2004), pp. 553-575. 

[2] M. AiNSWORTH AND H. A. Wajid, Dispersive and dissipative behavior of the spectral element 
method, SIAM J. Numer. Anal., 47 (2009), pp. 3910-3937. 

[3] D. Arnold, An interior penalty finite element method with discontinuous elements, SIAM J. Numer. 
Anal., 19 (1982), pp. 742-760. 

[4] A.K. Aziz and R.B. Kellogg, A scattering problem for the Helmholtz equation, in Advances in 
Computer Methods for Partial Differential Equations-Ill, vol. 1, 1979, pp. 93-95. 

[5] I.M. Babuska and S.A. Sauter, Is the pollution effect of the FEM avoidable for the Helmholtz 
equation considering high wave numbers?, SIAM review, 42 (2000), pp. 451-484. 

[6] I. Babuska, F. Ihlenburg, E.T. Paik, and S.A. Sauter, A generalized finite element method 
for solving the Helmholtz equation in two dimensions with minimal pollution. Computer methods in 
applied mechanics and engineering, 128 (1995), pp. 325-359. 

[7] I. Babuska and M. Suri., The h-p version of the finite element method with quasiuniform meshes.. 
Math. Modeling Numer. Anal., 21 (1987), pp. 199-238. 

[8] I. Babuska and M. Zlamal, Nonconforming elements in the finite element method with penalty, 
SIAM Journal on Numerical Analysis, 10 (1973), pp. pp. 863-875. 

[9] G.A. Baker, Finite clement methods for elliptic equations using nonconforming elements. Math. 
Comp., 31 (1977), pp. 44-59. 

[10] G. Bao and H. Wu, Convergence analysis of the pml problems for time-harmonic maxwell's equa- 
tions, SIAM J. Numer. Anal., 43 (2005), pp. 2121-2143. 

[11] J. -P. Berenger, a perfectly matched layer for the absorption of electromagnetic waves, J. Comput. 
Phys., 114 (1994), pp. 185-200. 



Pre- asymptotic Analysis of CIP-FEM and FEM 



24 



S.C. Brenner and L.R. Scott, The mathematical theory of finite element methods, Springer- 
Verlag, third ed., 2008. 

E. BuRMAN, A unified analysis for conforming and nonconforming stabilized finite element methods 
using interior penalty, SIAM journal on numerical analysis, 43 (2005), pp. 2012-2033. 

E. BuRMAN AND A. Ern, Stabilized Galerkin approximation of convection-diffusion-reaction equa- 
tions: discrete maximum principle and convergence, Mathematics of computation, 74 (2005), p. 1637. 

, Continuous interior penalty hp- finite element methods for advection and advection- diffusion 



equations. Math. Comp., 259 (2007), pp. 1119-1140. 

E. Burman, M.A. Fernandez, and P. Hansbo, Continuous interior penalty finite element 
method for Oseen's equations, SIAM journal on numerical analysis, 44 (2006), pp. 1248-1274. 

E. Burman and P. Hansbo, Edge stabilization for Galerkin approximations of convection- 
diffusion-reaction problems. Computer methods in applied mechanics and engineering, 193 (2004), 
pp. 1437-1453. 

Z. Chen and H. Wu, An adaptive finite element method with perfectly matched absorbing layers 
for the wave scattering by periodic structures, SIAM J. Numer. Anal., 41 (2003), pp. 799-826. 

P. G. CiARLET, The Finite Element Method for Elliptic Problems, North-Holland, Amsterdam, 
1978. 

P. CUMMINGS AND X. Feng, Sharp regularity coefficient estimates for complex-valued acoustic 
and elastic Helmholtz equations. Mathematical Models and Methods in Applied Sciences, 16 (2006), 
pp. 139-160. 

A. Deraemaeker, I. Babuska, AND P. BouiLLARD, Dispersion and pollution of the FEM solution 
for the Helmholtz equation in one, two and three dimensions. Int. J. Numer. Meth. Engng, 46 (1999), 
pp. 471-499. 

J. Douglas Jr and T. Dupont, Interior Penalty Procedures for Elliptic and Parabolic Galerkin 
methods. Lecture Notes in Phys. 58, Springer- Verlag, Berlin, 1976. 

J. Douglas Jr, JE Santos, and D. Sheen, Approximation of scalar waves in the space-frequency 
domain. Math. Mod. Meth. Appl. Sci, 4 (1994), pp. 509-531. 

B. EngquiST and a. Majda, Radiation boundary conditions for acoustic and elastic wave calcu- 
lations. Communications on Pure and Applied Mathematics, 32 (1979), pp. 313-357. 

B. EngquiST and L. Ying, Sweeping preconditioner for the helmholtz equation: moving perfectly 
matched layers, Multiscale Modeling and Simulation, 9 (2011), p. 686. 

X. Feng and H. Wu, Discontinuous Galerkin methods for the Helmholtz equation with large 
wave numbers., SIAM J. Numer. Anal., 47 (2009), pp. 2872-2896, also downloadable at 
|http : //arXiv . org/abs/0810 . 1475, 

, hp- discontinuous Galerkin methods for the Helmholtz equation with large wave number. Math. 

Comp., (2011, posted online). 

X. Feng and Y. Xing, Absolutely stable local discontinuous Galerkin methods for the Helmholtz 
equation with large wave number, Arxiv preprint arXiv: 1010.45631 (2010). 

B. Guo, Approximation theory for the p-version of the finite element method in three dimensions, 
part 1: approximabilities of singular functions in the framework of the jacobi-weighted Besov and 
Sobolev spaces., SIAM J. Numer. Anal., 44 (2006), pp. 246-269. 

B. Quo and W. Sun, The optimal convergence of the h-p version of the finite element method with 
quasi-uniform meshes., SIAM J. Numer. Anal., 45 (2007), pp. 698-730. 



Pre-asymptotic Analysis of CIP-FEM and FEM 



25 



[31] I. Harari, Reducing spurious dispersion, anisotropy and reflection in finite element analysis of time- 
harmonic acoustics, Computer methods in applied mechanics and engineering, 140 (1997), pp. 39-58. 

[32] U. Hetmaniuk, Stability estimates for a class of Helmholtz problems, Commun. Math. Sci, 5 (2007), 
pp. 665-678. 

[33] R.H.W HOPPE and B. WOHLMUTH, Element-oriented and edge-oriented local error estimators 
for nonconforming finite element methods, Math. Model. Numer. Anal, 30 (1996), pp. 237-263. 

[34] F. Ihlenburg and I. Babuska, Finite element solution of the Helmholtz equation with high wave 
number. I. The h-version of the FEM, Comput. Math. Appl., 30 (1995), pp. 9-37. 

[35] , Finite element solution of the Helmholtz equation with high wave number. H. The h-p version 

of the FEM, SIAM J. Numer. Anal., 34 (1997), pp. 315-358. 

[36] O.A. Karakashian and F. Pascal, A posteriori error estimates for a discontinuous galerkin 
approximation of second-order elliptic problems, SIAM Journal on Numerical Analysis, (2004), 
pp. 2374-2399. 

[37] J. Kim and D. Sheen, A priori estimates for elliptic boundary value problems with nonlinear 
boundary conditions, to appear, (1995). 

[38] J.M. Melenk, On Generalized Finite Element Methods, PhD thesis. University of Maryland at 
College Park, 1995. 

[39] J. M. Melenk and S. Sauter, Convergence analysis for finite element discretizations of the 
Helmholtz equation with Dirichlet-to- Neumann boundary conditions, Math. Comp., 79 (2010), 
pp. 1871-1914. 

[40] , Wavenumber explicit convergence analysis for Galerkin discretizations of the Helmholtz equa- 
tion, SIAM J. Numer. Anal., 49 (2011), pp. 1210-1243. 

[41] P. Oswald, On a BPX-preconditioner for pi elements. Computing, 51 (1993), pp. 125-133. 

[42] A.H. ScHATZ, An observation concerning Ritz-Galerkin methods with indefinite bilinear forms. 
Math. Comp., 28 (1974), pp. 959-962. 

[43] C. Schwab, P- and hp-Finite Element Methods, Oxford University Press, 1998. 

[44] J. Shen and L.L. Wang, Analysis of a spectral- Galerkin approximation to the Helmholtz equation 
in exterior domains, SIAM J. Numer. Anal, 45 (2007), pp. 1954-1978. 

[45] L.L. Thompson, A review of finite-element methods for time-harmonic acoustics, J. Acoust. Soc. 
Am., 119 (2006), pp. 1315 1330. 

[46] L.L. Thompson and P.M. Pinsky, Complex wavenumber fourier analysis of the p-version finite 
element method, Computational Mechanics, 13 (1994), pp. 255-275. 

[47] H. Wu, Pre-asymptotic error analysis of CIP-FEM and FEM for Helmholtz equation with high wave 
number. Part I: Linear version, to appear. 



